Unusual magnetic properties of the low-dimensional quantum magnet Na 2 V 3 07 
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We report the results of low-temperature measurements of the specific heat C P (T), ac suscepti- 
bility Xac(T) and 23 Na nuclear magnetic resonance NMR of Na2Va07. At liquid He temperatures 
C P (T)/T exhibits broad field- dependent maxima centered at To(H), which shift to higher temper- 
atures upon increasing the applied magnetic field H. Below 1.5 K the ac magnetic susceptibility 
Xoc(r) follows a Curie- Weiss law and exhibits a cusp at 0.086 mK which indicates a phase transi- 
tion at very low temperatures. These results support the previous conjecture that Na2Va07 is close 
to a quantum critical point (QCP) at fioH ~ T. The entire data set, including results of mea- 
surements of the NMR spin-lattice relaxation T 1 ~ 1 (T), reveals a complex magnetic behavior at low 
temperatures. We argue that it is due to a distribution of singlet-triplet energy gaps of dimerized 
V moments. The dimerization process evolves over a rather broad temperature range around and 
below 100 K. At the lowest temperatures the magnetic properties are dominated by the response of 
only a minor fraction of the V moments. 

PACS numbers: 75.40.-s, 75.30.kz, 76.60.-k, 75.30.Et, 
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I. INTRODUCTION 



Low-dimensional quantum magnets have recently been 
in the focus of both experimental and theoretical 
research 0, Q- Among different relevant materials of 
interest, oxides played a major role, in particular a vari- 
ety of vanadium oxides. A few years ago another mem- 
ber, Na2V3C>7, was added to this family by P. Millet 
and coworkers. The arrangement of the V magnetic mo- 
ments in this compound exhibits very peculiar geometri- 
cal features^ and in some respects, this material may be 
regarded as a spin S — 1/2 ladder- type compound with 
periodicity along the rung direction. This periodicity af- 
fects the topology of the arrangement of the magnetic 
moments. From the theoretical point of view it is still 
unclear to what extent this affects the magnetic proper- 
ties of the system |(j . 

Recent work has revealed that N^VaCv is an insulator 
with unusual low temperature properties pllll- Consider- 
ing the temperature dependence of the dc-magnetic sus- 
ceptibility x(T), it was arguedplQ that at temperatures 
T > 100 K, all the V ions adopt a tetravalent configura- 
tion (S = 1/2). The resulting moments act paramagnet- 
ically but are exposed to substantial antifcrromagnetic 
interactions. If we denote the ratio between an effective 
number of uncompensated V magnetic moments, each of 
spin S = 1/2, and the number of all the V ions as n mag , 
it follows that for T > 1Q0K, n mag — 1, and the para- 
magnetic Curie temperature Q p is of the order of -200 
K. Below 100 K, the V moments appear to be gradually 
quenched, such that below 20 K, only one out of nine V 
moments, exposed to a much weakened mutual coupling, 
contributes to x(T). In this regime, n mag = 1/9. From 
the 23 Na-NMR response of Na 2 V 3 07 it was concluded^ 



that in the presence of modest external fields H, up to 7 
T, and at temperatures of the order of 2.5 K and below, 
phase transitions occur at field dependent critical tem- 
peratures T a . Since T a shifts towards zero upon reducing 
H, it was suggested that Na2Va07 is close to a quantum 
critical point (QCP) at fi H «0T. 

Na2Vs07 crystallizes in a trigonal structure with a = 
10.886 A and c = 9.5380 A0. The space group is P31c 
and the unit cell contains 6 formula units. The V ions oc- 
cupy the geometric centers of O5 square pyramids, which 
in turn form the walls of nanotube-type structural ele- 
ments with a diameter of 5 A. These tubes are aligned 
along the c— axis of the crystal structure and the V ions 
occupy 9 sites on individual rings on the circumference 
of the tubes. The interactions among the V moments 
vary over a wide ranged . Generally, they involve some 
hybridization between V-neighboring O orbitals. The Na 
ions occupy four inequivalent sites. One of them, Nal, 
is located in the center of the nanotubes, and the rest, 
Na2, Na3, and Na4, is arranged around them. All the Na 
ions are bonded to oxygen ions which are at the vertices 
of complicated polyhedra around the Na ions. The same 
oxygen ions also form O5 square pyramids around the V 
ions. 

We report on measurements of the specific heat C P (T), 
ac-magnetic susceptibility Xac{T) and 23 Na-NMR on 
Na 2 V 3 7 . Our low-temperature C P (T)/T and C P (T) 
data reveal broad maxima at field dependent tempera- 
tures. Below 1.5 K, Xac{T) increases in a Curie-Weiss 
type fashion with decreasing temperature. A narrow 
maximum at 0.086K indicates a phase transition. Our 
new results support our previous claim that Na2Va07 
is close to a quantum critical point QCP at H = T. 
Anomalies in the temperature and field dependencies of 
the NMR response and x{T) at low temperatures indi- 
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cate different regions in the H — T phase diagram with 
different magnetic properties. It is argued that this re- 
flects a complex distribution of singlet-triplet energy gaps 
of dimcrizcd V moments. 



II. EXPERIMENTAL DETAILS 
A. Measurements 

The C p (T) measurements were made using a 
relaxation-type method at temperatures between 0.3 and 
40 K in either a 3 He refrigerator, or a 4 He-flow cryo- 
stat, both fitting into the bore of a 7 T superconducting 
magnet. For these measurements the sample consisted of 
small crystals, densely packed into a small thin-walled Cu 
container. In order to gain access to C P (T) of Na2Va07 
alone, the specific heat of the empty Cu container was 
measured separately. For the measurements of the ac- 
susceptibility between 0.018 and 1.8 K, we used an ac- 
impedance bridge with an rf-SQUID (superconducting 
quantum interference device) as the null detector. In 
this custom-built arrangement, the detection coils and 
the sample are placed inside the mixing chamber of a 
dilution refrigerator, in direct contact with the liquid he- 
lium mixture. For our measurements the sample was 
kept stationary inside the gradiometer-type arrangement 
of the detection coils [13. The ac driving field had an 
amplitude of 1.3 mOe and a frequency of 80 s _1 . 

The 23 Na NMR measurements between 0.1 and 300 K, 
were performed using standard spin-echo techniques and 
home-built spectrometers and probes. An unsuccessful 
attempt was made to observe the 51 V NMR signal at 
high (295 K) and low (near 1 K) temperatures. We con- 
cluded that at high temperatures, the spin-spin relax- 
ation time Ti for the 51 V nuclei must be very short, less 
than 5 [is. At low temperatures, the 51 V NMR signal 
may be distributed over a very broad range of frequen- 
cies. The central Zeeman transition (+1/2 <-» —1/2) of 
the 23 Na NMR spectra was obtained by integrating the 
spin-echo signal after a 7r/2 — tt, r /-pulse sequence, at 
a fixed applied magnetic field by changing stepwise the 
frequency and also at a fixed frequency by changing step- 
wise the magnetic field. The quadrupolar wings could not 
be observed except, perhaps, at very low-temperatures, 
where a weak NMR signal intensity was recorded around 
the central transition. Indeed, the results of a calcula- 
tion of electric field-gradients at the Na sites, using a 
point-charge model, confirm that for our randomly ori- 
ented powder-sample, the NMR signal intensity from the 
quadupolar wings is distributed over a broad frequency 
range. Therefore, all the NMR results presented below 
refer only to the central transition of the 23 Na spectrum. 
The spin-lattice relaxation time T\ was measured by first 
destroying the nuclear magnetization upon the applica- 
tion of a long comb of rf pulses, waiting for a variable 
recovery time t, and then measuring the integrated spin- 
echo intensity after a n/2 — n, r/-pulse sequence. For 
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FIG. 1: Cp/T as a function of T 2 for Na 2 V 3 7 at zero applied 
magnetic field. The solid line represents the sum expressed 
in Eq. 1 with the parameters quoted in the text. The broken 
line represents the lattice contribution. 

the NMR measurements below 20 K the sample was in 
contact with the 3 He- 4 He mixture of a dilution refriger- 
ator. At temperatures above 20 K, a 4 He flow cryostat 
was used to cover this temperature regime. Part of the 
NMR data were published previously in Ref. 8. 

B. Samples 

The material was prepared by growing very small crys- 
tals from melts of the starting composition Nai.gV2O5 
at reduced pressure. Densely packed material was used 
for measurements of C P (T) and Xac(T). For the NMR 
experiments the sample consisted of loosely packed pow- 
der. 



III. EXPERIMENTAL RESULTS AND 
RELATED CALCULATIONS 

A. Specific heat 

In Fig. 1 we display C p /T as a function of T 2 , mea- 
sured in zero applied magnetic field. Above 15 K the C p 
data is well represented by 

C'p^jT + Cde, (1) 

with 7 = 35 mJ/(mol K 2 ) and Cde, emphasized by the 
dotted line, representing the specific heat of the lattice 
due to acoustic phonons (Debye model with 0£> = 125 
K) and the lowest optical mode (Einstein model) at 88 
cm~ 1 [Tll|. The temperature variation according to Eq. 1 
is represented by the solid line in Fig. 1. At temperatures 
above 15 K, the specific heat is practically independent of 
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FIG. 2: C m /T = {Cp - Cde)/T as a function of InT for 
Na2V3C>7 measured in different external magnetic fields. The 
solid line emphasizes the logarithmic temperature variation. 



the magnetic field H, but at lower temperatures it varies 
significantly with increasing H . 

Considering the specific heat C m (T) = C P (T) — Cde, 
i.e., in excess of the usual lattice contribution of insula- 
tors, the parameter 7 turns out to be rather large. Since 
Na2V3C>7 is an insulator, the corresponding contribution 
to C m (T) is definitely not of (itinerant) electronic origin 
but most likely due to magnetic degrees of freedom. In 
Figs. 2 and 3 we display the excess specific heat C m , 
as defined above, plotted as C m /T versus T, including 
data for external magnetic fields ranging from to 7 
T. The temperature- and field-independent contribution 
registered above T = 15 K corresponds to the large 7T 
term in C P (T) mentioned above. 

In an applied magnetic field of 0.5 T we observe a loga- 
rithmic divergence at very low temperatures, emphasized 
by the solid line in Fig. 2, which represents 



C m {T)/T = -0.38 -ln(T/T K ), 



(2) 



where C m is measured in units of J mol _1 K -2 and 
Tk = 3.2 K. This behavior is identified between 0.37 K, 
the lower end of the covered temperature range, and 0.83 
K. A similar dependence, with the same pre-factor but a 
lower Tk characterizes the C m (T)/T data measured in 
zero field. Logarithmic divergences of C m /T at low tem- 
peratures were often reported for a variety of fermionic 
systems close to QCPsjll E E3- 

The C m (T) data, obtained for fields exceeding 1 T, 
however, exhibit field dependent maxima in C m (T)/T 
and, even for H — T, in C m (T), as demonstrated in 
Figs. 2, 3 and 4. The temperatures of the maximum for 
C m (T) does not vary much upon changing the magnetic 
field H, at least for fioH < 6 T. The position of the maxi- 
mum of C m (T) /T, however, shifts to higher temperatures 
with increasing field (see Fig. 5). The solid lines in Figs. 
3 and 4 represent the best fits to the data using a model 



FIG. 3: Cm/T as a function of T for Na2Vs07 measured in 
different external magnetic fields. The solid lines represent 
the best fits to the data using a model described in the text. 
The horizontal broken lines represent the constant 7 term (see 
Eq. 1). 




FIG. 4: C m as a function of T for Na2Vs07 measured in 
different external magnetic fields. The solid lines represent 
the best fits to the data using a model described in the text. 
The broken lines represent the 7T term (see Eq. 1). 



to be described below. The model is chosen such that the 
calculation of C m (T) and C m {T)/T provides an accept- 
able agreement with experiment. The data below 20 K 
(see Figs. 1, 3 and 4) suggest that C m (T) = jT+C st (T), 
whereby C s t(T) contributes significantly only below 15 
K. In the following we concentrate on C s t{T) and the 
7T term will only be considered again in the discussion 
section. 

It turned out that the best agreement between experi- 
ment and model calculation could be obtained by assum- 
ing that C s t(T) arises from excitations in a collection 
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FIG. 5: T m ax, the position of the maxima of (C m (T) — jT) 
(filled circles) and C m {T)/T (filled triangles), as a function 
of H for Na2Vs07. The circle with the center dot represents 
the sharp peak in Xac(T) (see Fig. 8) and the empty circles 
the position of the maxima at T a in the NMR spin-lattice 
relaxation rate. 



FIG. 6: Probability distribution P(A) = Pi (A) + P 2 (A) of 
the spin singlet-triplet energy gap A. P(A) was obtained 
as the sum of two Gaussians from fitting the C m (T) and 
C m (T)/T data. Inset: Magnetic field dependence of the cen- 
ters (T m axi and T max 2) and of the widths (Wi and W2) of 
the Gaussians (Pi and P2). 



of dimerized pairs of V moments, i.e., an ensemble of 
singlet-triplet units with a distribution P(A) of singlet- 
triplet energy gaps A. A more detailed account on the 
justification of this choice is given in the discussion sec- 
tion. We assume that P(A) is the sum of two Gaussians, 
Pi (A) and P 2 (A). Pi (A) is narrow with a width Wi/k B 
of the order of 4 K [HWHM = 0.83 x W 1 /k B = 3.2 K) 
centered at T max i = 0K and P 2 (A) is a broader Gaus- 
sian centered at T maX 2 ~ 9 K with a width W^/fce of 
about 10 K {HWHM = 0.83 x W 2 /k B = 8.3 K). The 
distribution of P(A) is thus given by 



P(A) 
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exp • 



exp • 
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W 2 



}■ (3) 



Keeping in mind that a formula unit of the substance 
contains 3 V ions or 3/2 V pairs, and that at low temper- 
atures, n rnag = 1/9 as explained in the introduction, the 
prefactors are chosen to fix the total population (area) 
to (l/9)(3/2) = 1/6 pairs per formula unit. The ratio 
between the relative populations of the narrow and the 
broad Gaussian is about 0.5 and is the same for all the 
data sets. Taking W%, T maX 2 and W2 as free param- 
eters (T max i is fixed to zero) varying slightly with H, 
we can reproduce very well all the available C m (T) and 
C m (T)/T data as demonstrated by the solid lines in Figs. 
3 and 4. For this we have calculated the specific heat due 
to a collection of individual singlet-triplet units in an ap- 
plied magnetic field H. The levels are considered to be 
(5-peaks with the energy of the singlet level at zero and 



of the triplet levels at A, A + 2[IbH and A — 2/i#iJ, re- 
spectively. The calculations are based on the distribution 
function P(A) given in Eq. 3. 

P(A) in zero magnetic field, given by Eq. 3 and result- 
ing from inserting the best fit parameters, is represented 
in Fig. 6. The inset displays the magnetic- field induced 
variations of the parameters W± , T max 2 and W2 ■ Ideally 
they should be field independent, but the best agreement 
with experiment is obtained by invoking a weak field de- 
pendence which, most likely, reflects the limits of our 
model. Inspection of Fig. 6 reveals a small fraction of 
singlet-triplet units with negative gaps, A < 0, implying 
that for a minority of pairs the spin-triplet configuration 
is energetically more favorable. This would be rather sur- 
prising if the inter-ring exchange was a simple superex- 
change between appropriate d xy orbitals, as assumed 
for instance in Ref. |15| . because the couplings would 
then all be antiferromagnetic. But since Mazurenko and 
coworkers)^, on the basis of an LDA+U estimate of the 
exchange integrals, conclude that inter-ring couplings are 
ferromagnetic, the quoted implication is actually plausi- 
ble. The model accounts for the field-dependent maxima 
of both C P (T) and C P (T)/T and for the substantial in- 
crease of C P (T)/T with decreasing temperatures in ex- 
ternal magnetic fields of and 0.5 T, respectively (see 
Figs. 2 and 3). 

The above discussion implies that the anomalies of 
C P {T) and C P (T)/T do not reflect cooperative phase 
transitions. However, the low-temperature anomalies at 
T a in the NMR spin-lattice relaxation T 1 _1 (T), reported 
in Ref. 8, reveal changes in the dynamics of the V mo- 
ments in the same temperature range. We address this 
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FIG. 7: Normalised magnetic entropy (S m — yT) per mol 
of V ions as a function of temperature for different external 
magnetic fields. The horizontal broken lines represents 1/9 
of the normalized total entropy, Stotai/ (3 R In 2) = 1. Inset: 
Temperature dependence of the magnetic entropy including 
the contribution of the 7T term for [10H = 3T. 

issue in the discussion section. Since the maxima in 
C p (T)/T or in C P (T) are well accounted for by our model 
and no other specific-heat anomalies are discernible, also 
the features in T 1 _1 (T) at T a cannot reflect transitions to 
conventional magnetically-ordered phases. 

The magnetic entropy S m (T) was obtained by numeri- 
cal integration of C m (T) /T and asserting that S m (0) = 0. 
In Fig. 7 we display S m (T) — 7T for magnetic fields of 0, 
0.5, 3, 5 and 7 T. These data imply that only a very small 
magnetic entropy S m (T) — 7T is released at low temper- 
atures, reaching approximately 1/9 of 3i? • In 2 at 20 K 
or above. This is the value expected for the entropy of a 
system of dimers that is described by Eq. 3, and it im- 
plies that indeed, n mag = 1/9 below 20 K. This observa- 
tion adds convincing support for our previous conjecture 
that the V magnetic moments in Na2V3 07 are gradu- 
ally quenched with decreasing temperature between 100 
and 20 K. As may be seen in the inset of Fig. 7, above 
9 K, S m (T) increases almost linearly with T, reflecting 
the importance of the large 7T term of the specific heat. 
As mentioned before, the significance of this term is ad- 
dressed in the discussion section. 



B. ac-susceptibility 

In Fig. 8 we display the results for Xac(T), after sub- 
traction of a T— independent signal. Between 0.11 and 
1.5 K the data is well represented by 

Xac(T) = X o + C/(T-e p ). (4) 

The very small paramagnetic Curie temperature P = 
—0.084 K reveals a very weak coupling of the remaining 



FIG. 8: Temperature dependence of the ac-magnetic suscep- 
tibility of Na2Vs07 in zero applied magnetic field. The solid 
line represents the best fit to a Curie- Weiss law. The sharp 
peak in Xac(T) at 0.086 K reveals a phase transition at an 
unexpectedly low temperature. Inset: Xac vs - T. 

magnetic moments of V (see also the inset of Fig. 8). 

Next we estimate the concentration of the V magnetic 
moments that contribute to x{T). Since the absolute 
value of Xac cannot be extracted directly from the Xac 
measurements, we simply match the Curie- Weiss parts 
of Xac(T) and X dc (see Ref. 8) at T = 2 K. From this 
procedure we find that the low-temperature Curie con- 
stant C of Eq. 4 is only half of the corresponding value at 
temperatures above 2 K. This is interpreted as reflecting 
a further reduction of n mag by the same factor. In sum- 
mary, from all the results for x(^) we conclude that in 
Na2V3C>7 near room temperature, all the V moments are 
active and n mag — 1. Between 2 and 20 K, n rnag = 1/9 
and between 0.11 and 1.5 K, n mag = 1/18. 

The very sharp peak in Xac(T) at T max (H w 0) = 0.086 
K signals a phase transition involving spin degrees of free- 
dom. Other than a slight suppression of the amplitude of 
the peak, no field- induced changes are observed in Xac{T) 
measured in small fields of the order of 150 Oe (data not 
shown). Since at these temperatures the concentration 
of active magnetic moments is very small (n m ag = 1/18), 
it is again very unlikely that this transition is of conven- 
tional ferro- or antiferromagnetic type. Most probably it 
reflects the freezing of the few remaining active magnetic 
moments of the V ions or, more accurately, it reflects the 
freezing of dimerized V ions with a very small energy gap 
or with a triplet ground state. 



C. NMR at low temperatures, T < 40 K 

In Fig. 9 we display examples of the 23 Na-NMR cen- 
tral transition signal at very low temperatures for three 
fixed frequencies vl of 11.03, 31.73 and 77.5 MHz, corre- 
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FIG. 9: Examples of the 23 Na-NMR central line of Na 2 V3O7 
for three different frequencies of 77.5, 31.73 and 11.02 MHz 
at very low temperatures. 




T7v L (K/MHz) 

FIG. 10: 23 Na-NMR linewidth F s of Na 2 V 3 07 as a func- 
tion of T/pl below 40 K, with ul denoting the resonance 
frequency. The vertical dotted line at T a /vL separates the 
temperature independent- from a strongly temperature de- 
pendent regime. The solid line represents a Curie-type law 
and the arrows mark the temperatures where T 3 (T) departs 
from the high-T behavior. 



sponding to resonance fields fi H of 0.98, 2.825 and 6.9 
T, respectively. The comparison of the central transition 
signals for vl = 31.73 MHz at 17 and 0.1 K emphasizes 
the very large increase of the linewidth FWHM = T s with 
decreasing temperatures below 20 K (see also Ref. 8). 
Since in this temperature regime, xC^O follows a Curie 
type behavior with a small paramagnetic Curie temper- 
ature Q p « -2 K, the high temperature (T > 8 P ) NMR 
linewidth for a randomly oriented powder is expected to 
vary as T s oc M oc H/T. In order to check this pre- 
diction we plot the linewidth of the 23 Na-NMR central 
line of Na2Va07 in Fig. 10 for three different frequen- 
cies of 77.5, 31.73 and 11.02 MHz as a function oiT/v L , 
with vi, denoting the corresponding resonance frequency 
(vl oc H). We note that, below 40 K, the initial in- 
creases of the linewidth with decreasing temperature fall 
onto the same curve for all data points, indicating that 
indeed T s oc (T / H)^ 1 . This is emphasized by the solid 
line which represents a Curie type temperature depen- 
dence. This behavior breaks down at low temperatures 
and in Fig. 10 the onset of deviation is emphasized by 
the vertical arrows. Considering that T s measures the 
width of the distribution of local susceptibilities Ax, our 
data suggest that Ax is large at low temperatures. 

In Fig. 10, another feature worth to be noticed 
emerges. Below a constant value T/i>l « 0.02 (K/MHz), 
the linewidths T s are temperature independent, at a level 
that changes significantly with applied magnetic field. 
The vertical solid line at Tjv^ = T a /vL represents ap- 
proximately the positions of the maxima in T 1 ~ 1 (T) at 
T a (H) reported in Ref. 8. As we already mentioned 
above, it is not justified to introduce a phase boundary 
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FIG. 11: Temperature dependence of the NMR spin-lattice 
relaxation rate in applied magnetic fields of 2.82 and 6.81 T 
at temperatures below 3 K. Prominent anomalies define T a . 

at T a (H). As previously described in Ref. 8, Tf^T) 
exhibits yet another anomaly in the form of an abrupt 
change of slope at Tb(H), distinctly below T a (H). For 
the reader's benefit, we partly recall these features in 
Fig. 11. The slope change in T^ 1 was interpreted as 
indicating the formation of a small gap in the spectrum 
of V spin excitations at Tb(H). In view of lacking addi- 
tional evidence from specific-heat data, we refrain from 
a further discussion of this feature in the present work. 

Our T s data are difficult to reconcile with scenarios of 
either a simple type of field-induced magnetic ordering or 
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FIG. 12: 23 Na-NMR linewidth (T s ) of the Zeeman central 
transition of Na2Vs07 as a function of T. The data was 
measured at 66.25 MHz and the solid line emphasizes the 
Curie- Weiss behavior. Inset: NMR line shift if as a function 
of T for data measured at 38.7 MHz. 

a spin freezing at T a (H). In such cases, the linewidth at 
very low temperatures would be caused by internal static 
fields and hence would not depend on H, contrary to our 
observations. If, nevertheless, "spin freezing" is assumed 
to occur at T a then, only a minority of spins is involved 
in the phenomenon, as suggested by the results for the 
specific heat. These particular spins are also dimerized as 
all the others, but either the related singlet-triplet energy 
gaps are small or zero, or a triplet ground state is estab- 
lished. This may be plausible, but does not address the 
striking fact that T a varies linearly with H . Therefore, 
we abandon this scenario at the outset. 



D. NMR at high temperatures T > 40 K 

Next we examine the 23 Na-NMR data obtained in the 
temperature region above 40 K. They provide additional 
experimental evidence for the drastic change in the V 3d 
electron system around 100 K, as previously concluded 
from x{T) data in Ref. 8. The first evidence is pro- 
vided by the temperature dependence of T s as shown 
in Fig. 12. The Curie- Weiss type enhancement with de- 
creasing temperature, emphasized by the solid line in Fig. 
12, is abruptly terminated at approximately 100 K. The 
plateau region between 100 and 50 K is compatible with 
a loss of free ionic V moments and hence a reduction of 
n ma g- The onset of yet another rapid increase of T S (T), 
with decreasing T is consistent with a very much reduced 
Q p for temperatures below 40 K. 

We emphasize that at low temperatures, the width of 
the 23 Na-NMR line is very large and is clearly dominated 
by the magnetism of the V moments. On the other hand, 
the relative shift K of the 23 Na-NMR line is rather mod- 
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FIG. 13: Temperature dependence of the NMR spin-lattice 
relaxation rate in an applied magnetic field of 5.88 T at tem- 
peratures above 40 K. A prominent maximum defines Tab- 



est and the magnetism of the V moments may only be 
one of several factors affecting K(T). Nevertheless, we 
note a considerable enhancement of K(T) from T > 150 
K to T < 75 K (see inset of Fig. 12). 

In Fig. 13 we show T-^ (T) measured in an applied 
field of 5.88 T. Above 200 K, Tf 1 (T) is approximately 
T-independent, as expected for a relaxation driven by 
the flips of paramagnetic moments. At lower temper- 
atures, Tf 1 (T) passes through a prominent maximum 
at Tab near 100 K. Since the relaxation is mainly due 
to the flips of V moments 8], which as mentioned above 
should result in a T— independent relaxation, this max- 
imum of Tf ^T) must be associated with changes in the 
dynamics of the V moments. From specific-heat and 
magnetic-susceptibility data we infer that the most likely 
scenario involves the formation of V-V dimers of the V 
moments. In this case the loss of magnetic degrees of 
freedom associated with the formation of the dimers, 
i.e., with the reduction of n mag , may be thought of as 
a gradual reduction of the fluctuation frequencies of the 
involved moments. This process is expected to result in a 
broad peak in Tf (T) similar to observations in precur- 
sor regions of magnetic ordering and spin freezing. The 
recent suggestion [2(j that the reduction of n mag maybe 
explained by conventional crystal-field effects and spin- 
orbit coupling of the V 4+ ions is judged as unconvincing, 
because in such a scenario, the peak in Tf (T) is not 
expected. 



IV. DISCUSSION 

A. The different regions in the H — T plane 

In this section we summarize the obtained experimen- 
tal results for Na2V 3 C)7 and, as before, we characterize 
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different regions in the H~T plane by n mag , which repre- 
sents the concentration of those moments that dominate 
the increase of x with decreasing temperature. 

The results for the specific heat, the NMR response (in- 
cluding NMR data from Ref. 8) as well as the ac— and 
dc— susceptibility, imply several characteristic tempera- 
tures that are indicated by anomalies in the temperature 
dependences of 

1. T-j - , exhibiting maxima at Tab and T a , 

2. r s , revealing changes in the T dependence at Tab 
and at T a , 

3. C m (T)/T, and C m (T), with maxima at field- 
dependent temperatures, 

4. Xac, with a peak at T max (H = 0) and 

5. Xdc, implying two different temperature regimes 
with different n mag (see also Ref. 8). 

Well above 100 K, in a region that we denote as re- 
gion A, all the V magnetic moments contribute to the V- 
based 3d electron susceptibility, i.e., n mag = 1. The low- 
temperature boundary of this conventional paramagnetic 
phase is marked by a very prominent and broad peak in 
Tf 1 (T) at T A b ~ 100 K (see Fig. 13). Around the same 
temperature we note significant changes in the tempera- 
ture dependences of the NMR linewidth T S (T), the NMR 
line shift K(T) and x(T). Most of these features are con- 
sistent with a gradual reduction of n mag upon decreasing 
temperature. 

The observed variation of the chemical shift K across 
100 K (see inset of Fig. 12) most likely reflects the on- 
set of distortions of the oxygen polyhedra surrounding 
the Na ions at temperatures between 75 and 125 K. This 
feature and the concomitant reduction of n mag suggest 
that a dimerization of the V magnetic moments is the 
most likely process to account for the change of the mag- 
netic properties from region A to region B, the dimerized 
phase. At temperatures well below Tab, i-e., between 20 
and 2 K, the spin system is characterized by n mag = 1/9. 

In the same temperature regime, i.e., well within re- 
gion B, we note an excess specific heat and part of it 
varies linearly with T. Inevitably, it has to be related 
to excitations in the magnetic subsystem of this insu- 
lating compound. Recalling the evidence that "effec- 
tively" still 1/9 of all the V ions act paramagnetically 
in this regime, it is tempting to ascribe this 7T term 
to S = 1/2 spinon excitations that occur in antifcrro- 
magnetic Heisenberg chains in the low-temperature limit. 
This interpretation does not make sense, however. In this 
situation [l^, CjJ Us LUJ 

7 = C p /T = 2N A k%/(3J), (5) 

where T <C J/ks, and J is the intra-chain exchange inter- 
action. In view of the reduced value of n mag , the experi- 
mental value of 7 implies that J/fcs w9K, incompatible 
with the request that T <C J/ks for Eq. 5 to be valid. 



We emphasize once more that the anomalies in T[~ (T) 
at T a (see Fig. 5), the maxima of C m (T, H) and those of 
C m {T,H)/T do not reflect new phase boundaries. The 
situation is insofar complicated that in zero magnetic 
field and at temperatures below 1.5 K, region B is char- 
acterized by n mag = 1/18. Thus the cusp in x{T) at 
T m ax(H « 0) = 0.086 K most likely indicates a spin- 
freezing type phenomenon. In view of the very low con- 
centration of active V moments and, by consequence, the 
randomness of their mutual interaction, it seems rather 
unlikely that a state of conventional magnetic order is 
adopted. As a possible scenario we suggest that the freez- 
ing at T max only invokes spins that occupy the very few 
triplet ground states that are compatible with the Pi (A) 
distribution of our model. The anomalies in T 1 _1 (T) at 
T a most likely reflect temperature-induced changes in the 
dynamics of the V moments. Note that in this case the 
dimers of P2(A) play only a minor role, since the vast 
majority of them have singlet ground states with energy 
gaps relatively large compared with those of Pi (A). For 
applied fields of the order of 0.5 T or larger, however, 
the distribution of minimum excitation gaps changes con- 
siderably. The dimers (with a triplet ground state for 
HoH = T) acquire now a gap of 2/i^iJ and play no role 
here. This suggests yet another phase, different from 
phase B and stable only at very low temperatures and in 
very small or zero external magnetic fields. 



B. The distribution of V-V interactions 

This subsection is devoted to a discussion of our model 
assumptions, the calculation of the low temperature spe- 
cific heat and to an attempt to understand the complex 
magnetic behavior of Na2Va07 at temperatures close to 
T = K. 

First of all we explore the implications of the crucial 
observation that well below Tab only 1/9 of the V mo- 
ments contribute significantly to the Curie- Weiss type 
susceptibility and that the rather small paramagnetic 
Curie temperature 9 p w —2K implies only weak inter- 
actions between them. We start by assuming that 8/9 
of the localized V spins are strongly dimerized, i.e., with 
rather large singlet-triplet energy gaps, and hence with- 
out much influence on x(T) below 20 K. The rest, 1/9, 
is first considered to be free. In this scenario, C s t{T) at 
low temperatures would exclusively be due to the free 
magnetic moments and therefore, it should be well rep- 
resented by a contribution due to a collection of identical 
two-level systems, each with a field-dependent gap given 
by A = 2[j,bH, where H is the applied magnetic field. 
This simple model, however, cannot account for the ex- 
perimental data. The gap vanishes for H = T, and no 
peak in C st (T) in H = T is expected, vastly different 
from what we observe (see Fig. 4). Other simple mod- 
els which assume, for instance, that the spin system is 
formed by a collection of identical singlet-triplet units, re- 
placing the pairs of free moments and leaving the singlet- 
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triplet energy gap as a single- value free parameter, also 
fail to account for the observed C m (T) data. This situa- 
tion clearly suggests that the considered V spins (1/9 of 
the total) can neither be assumed to be free nor their local 
environments to be identical. This should be contrasted 
with the properties of the effective spin-chirality model 
of odd- leg spin tubes [2jJ in which the ground state is 
spontaneously dimerized and the effective S = 1/2 spins 
of the rings are frozen in a regular pattern j^. Finally 
we recall that also the 7T term in C m (T) needs to be 
accounted for. 

The previous conclusions from x(T) and the Schottky- 
type features of C s t (T) suggest that the spin system may 
be represented by a collection of dimerized pairs of V 
moments with a distribution P(A) of singlet-triplet en- 
ergy gaps A. Since the situation is rather complex we 
proceeded in two steps. In a first step we again did not 
consider the 7T term in C m (T), i.e., we neglected the 
contribution to the specific heat due to the above men- 
tioned dimers with large singlet-triplet gaps (8/9 of the 
total). The remaining relevant dimers were assumed to 
exhibit gaps with the distribution P(A) given in Eq. 3. 
The best fits to the data turned out to be consistent with 
rather modest values for the singlet-triplet energy gaps. 
Physically, a distribution P(A) of energy gaps is not un- 
reasonable because, according to Ref. 9, a large number 
of different V-V couplings and, in addition, frustration 
effects may be important to understand the magnetic 
features of Na2Va07. Our model is based on a com- 
plete dimerization of all V spins at low temperatures. 
Most of the dimers (8/9) are characterized by fairly large 
singlet-triplet energy gaps, such that their influence on 
x(T) at very low temperatures may be neglected at this 
stage. This forces the minority of spins, 1/9 of the to- 
tal, to adopt dimerization partners from different rings, 
the characteristic structural element mentioned in the 
introduction, via weaker interactions. Therefore, due 
to the geometrical constraints of the structure and the 
expected wide range of energy characterizing the cou- 
plings between the V moments, a distribution P(A) of 
singlet-triplet energy gaps is likely to occur as the re- 
maining V spins dimerize. In this sense we argue that 
the model described in section III. a (see Eq. 3) which is 
successful in reproducing the experimental C s t (T, H ) and 
C s t(T, H)/T data, may be justified in the way sketched 
above. Naturally the same model is expected to repro- 
duce x(T) and we calculate the magnetic susceptibility 
Xi,2(T) due to the dimerized spins with the distribution 
of singlet-triplet gaps given by Eq. 3. The result is dis- 
played in the main frame of Fig. 14 as the dotted line. 
The result of this calculation, although reflecting a Curie- 
Weiss type behavior of the susceptibility, is in substantial 
disagreement with the measured data, shown as open cir- 
cles in the same figure. This is not surprising, however, 
because we have, as pointed out above, so far neglected 
the contribution due to the strongly dimerized V spins. 

Therefore, in the next step we attempted to identify 
the origin of the 7T term of the specific heat and the 
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FIG. 14: Inverse susceptibility vs - T. The circles rep- 
resent the data taken from Ref. 8. The dotted line is the 
contribution X\ \ due to the distribution P(A) of spin singlet- 
triplet gaps (see text) represented in Fig. 6. The solid line 
is obtained by including the additional contribution to x(T), 
X3, from the distribution P3 (A) of singlet-triplet gaps (see 
text). The inset display the resulting contribution of Ps(A) 
to the specific heat. 



missing contribution to x(P)- Following our previous 
procedure, we assumed that also the impact of dimerized 
spin pairs with larger singlet-triplet gaps can be captured 
with a distribution of gaps, here denoted as p3(A). More 
or less by trial and error and seeking agreement with the 
experimental data, we established Pj(A) to be given by 

( 1.26 x 10" 3 , -15K < A/k B < 250K 
P 3 (A) = I 1.00 x 10~ 2 , 250K < A/k B < 350K 
[ 0, otherwise. 

(6) 

The prefactors have been adjusted, such that the to- 
tal population (integrated area) represents (8/9) x (3/2) 
dimerized pairs of V moments per formula unit of 
Na2V30y. Note that from Eq. 6 it follows that about 
3% of the gaps A are of the order of 20 K or less, i.e., 
are not "strongly dimerized". With this choice of Pj(A) 
the model yields both, the necessary correction to Xi2(P) 
in order to reproduce fairly well the experimental data 
for x(T) (see Fig. 14), as well as the contribution to the 
specific heat which, in good approximation, reproduces 
the observed 7T term (see the inset of Fig. 14). We em- 
phasize that we do not claim that the above described 
Ps(A), nor P(A) given in Eq. 3, are the only possible 
or even the optimal choices for characterizing the dimers. 
Nevertheless, this comparison between experimental data 
and model calculation gives strong support for the claim 
that at temperatures well below 100 K, all the V ions are 
dimerized. 

In summary, the complex magnetic properties of 
Na2V3C>7 may be understood by assuming that at tem- 
peratures well above 100 K the system responds param- 
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agnetically with the localized V magnetic moments in- 
teracting with each other on a wide range of coupling 
energies. At temperatures of the order of 100 K, the V 
moments form pairs with mostly singlet (A > 0), but 
also with some triplet (A < 0) ground states. The com- 
parison with experiment suggests that the singlet-triplet 
gaps A adopt values with a broad a distribution P(A). 
The total distribution P(A) contains at least three terms, 
which here are represented by two Gaussian functions, as 
described in Eq. 3, for 1/9 of the V moments, and an 
additional very broad distribution, described in Eq. 6, 
which characterizes the rest of the V moments (8/9) . The 
evidence that leads to the claim that the low-temperature 
magnetic features are dominated by only 1/9 of the mo- 
ments is most clearly provided by the respective entropy 
release shown in Fig. 7. 

The model calculation captures the variation of the 
magnetic properties over a broad temperature range 
around and below Tab- The peak in T 1 _1 (T) at Tab re- 
veals that the dimcrization process must be accompanied 
by dramatic changes in the dynamics of the V moments. 
Further changes in the dynamics of the V moments, this 
time involving a slowing down of the fluctuations between 
the occupied energy levels of the dimers, seem to occur 
at much lower temperatures, of the order of T a , as judged 
by the anomalies in T 1 _1 (T) shown in Fig. 11. In a first 
approximation we argue that in this regime, T 1 _1 (T) is 
basically determined by transitions between spin levels of 
individual dimers with low energy gaps. The relaxation 
due to a collection of identical V dimers is given bv|23| 



At 



(1 + L0 2 T 2 )' 



(7) 



with r as the correlation time of the dimer spin transi- 
tions, ui = 2ttvl with i/£ oc H as the Larmor frequency 
and A is a constant which includes the hyperfine-field 
coupling between the considered dimers and the probed 
nuclei. Assuming that the T^ 1 variation is predomi- 
nantly due to a variation of t, Eq. 7 provides a max- 
imum of T-f^r) for lot = 1. It also seems reasonable 
to assume that the transitions are thermally activated. 
In the simplest case, the transitions are activated over a 
single barrier U and 1/t oc cxp(— U/ksT). In case of a 
distribution f(U) of barriers U, an effective correlation 
time T e f f is given by 



1 



r eff 



oc / exp 



U 
k^T 



f(U)dU. 



(8) 



If f(U) is assumed to be constant for the relevant values 
of U, this reduces to 



1 



exp 



U 



dU ocT. 



(9) 



Eqs. 7 and 9 explain a strong increase of T x with 
decreasing temperature in the limit wT e ff << 1 where 
T^ 1 w At oc T -1 . Likewise the decrease of Tf 1 beyond 



the maximum is also consistent with our assumptions, 
because T^ rs A/{lo 2 t) oc T for Lor e ff >> 1. Since the 
maximum in Tj - (T) is tied to the condition lot = 1, it 
follows that T a = T max oc 1/t oc lo oc H, as observed ex- 
perimentally. We are well aware that Eq. 7 only crudely 
captures the characteristics of our spin system. Without 
more detailed information on the spin system, however, 
we cannot offer a more sophisticated description of the 
dynamics of the dimers. Nevertheless, the fact that the 
data is consistent with r e 1 1 oc 1 /T suggests that the spin 
dynamics is rapidly slowing down with decreasing tem- 
peratures. As may be seen in Fig. 5, T max w 0.2 K for 
[i H s» 1 T. This field value corresponds to a frequency 
of 11.02 MHzj3|, a rather low value. The extrapolation 
of T a (H) = T max (H) shown in Fig. 5 suggests again that 
the system is close to a QCP, as previously claimed -8j. 

Finally we point out that one aspect of our data, 
the temperature independent but field dependent NMR 
linewidth T s at low temperatures (see Fig. 10), is not 
explained by our model. Commonly T S (T,H) is pro- 
portional to the distribution of local fields A/i; oc at the 
probed nuclei. As demonstrated in Fig. 10, this is also 
true in our case at temperatures exceeding T a (H) and 
consistent with our model which implies a monotonic in- 
crease of the total magnetization M with decreasing tem- 
perature. In this sense, the abrupt change of behavior of 
T s at T a /i/L is unexpected, indicating that our model 
misses an important ingredient which is effective only at 
low temperatures. The relevant dimers (with small gaps) 
in Na2V3C>7 are expected to be rather extended in real 
space. The implicit model assumption of non interacting 
dimers may thus be a rather crude approximation. A 
more detailed knowledge of the couplings between the V 
moments seems to be needed for a self consistent inter- 
pretation all the experimental data of Na2Va07. 



V. CONCLUSION 

The magnetic properties of Na2Va07 are qualitatively 
different in different regions in the H — T plane. These 
differences may be thought of as the result of a gradual 
"quenching" of the localized V moments with decreas- 
ing temperatures below approximately 100 K. But from 
the analysis of the specific heat we conclude that this 
" quenching" reflects a process of dimerization at temper- 
atures of the order of 100 K, resulting in a distribution 
P(A) of singlet-triplet energy gaps A. For ID systems, a 
simple dimerization generally leads to observable changes 
in the crystal structure. Such changes have been searched 
for, but were not observed in optical studies [Tl|. From 
our results it is suggested that the V — V couplings in- 
volve a broad range of interactions and given the unusual 
structural features of Na2Vs07, which may include some 
randomness, it is not clear that changes in the structure 
that could be detected in scattering and/or optical exper- 
iments should be expected here. We note, however, that 
the temperature induced variation of the 23 Na chemical 
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shift shown in the inset of Fig. 12, may well be due to 
subtle structural changes accompanying the compensa- 
tion of the moments. 

Anomalies in T^ l {T,H) and a cusp in x{T) at 0.086 
K reveal that Na2Va07 is close to a quantum critical 
point, as previously claimed. Finally the broad distribu- 
tion of couplings suggests that the relevant low-energy 
effective magnetic model should contain some level of 
randomness |3 HH Ull- Given the complexity of the 
V network, it is at present impossible to say whether 
these properties, including the behavior of T S (T,H) at 
low temperatures, may consistently be explained with a 
purely magnetic model with some coupling to the lattice, 
or whether some extrinsic disorder, for instance due to 



some missing Na +1 ions in the structure, has to be taken 
into account. 
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